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Abstract 



We present a measurement of the slope parameter a for the 77 — > 3tt° decay, 
with the KLOE experiment at the DA$NE 0-factory, based on a background 
free sample of ~ 17 millions rj mesons produced in radiative decays. By 
fitting the event density in the Dalitz plot we determine a = —0.0301 ± 
0.0035 stat iaoo35 syst . The result is in agreement with recent measurements 
from hadro- and photo-production experiments. 

Keywords: e + e~ collisions, 4> radiative decays, r\ decays 



1. Introduction 

The decay r] — > 3tt, 7r + 7r~7r°and 37r° , though is a major decay mode 
of the r) meson, violates isospin symmetry. Since contributions from the 
electromagnetic interaction are strongly suppressed by chiral symmetry [l[ 
this decay is mainly due to the isospin breaking part of the QCD Lagrangian: 



so that in principle it offers a way to determine the mass difference of the 
up-down quarks. Moreover, the selection rule AI = 1 allows us to relate the 
amplitudes for the two decays using isospin symmetry: 



Theoretical predictions for the decay amplitude have been obtained in the 
framework of Chiral Perturbation Theory (ChPT): the low energy effective 
field theory for QCD. Leading order (LO) ChPT predictions [2| based on 
current algebra underestimate the i] decay rates by a factor of ~ 4. One 
loop (NLO) calculations which include the n — ir rescattering effects 0] pre- 
dict higher rates but still below the observed values. Some improvements 
are obtained by computing unitary corrections [4] to NLO usinga dispersion 
relation for the decay amplitude derived by Khuri and Treiman |5j. Recently, 
more advanced calculations have become available. In reference j6[ the au- 
thors use U(3) ChPT, in combination with a coupled channels method, and 
treat final state interactions by means of the Bethe Salpeter equation ob- 
taining good agreement with measured decay widths and spectral shapes. In 
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Reference [7J a full NNLO computation is performed showing sizable correc- 
tions to the NLO result. 

The Dalitz plot of a three body decay is described by two kinematical vari- 
ables which for three identical particles in the final state, reduce to a single. 
In the rj — > 37T° decay this variable is chosen by convention to be: 

z= 2_ ' t3E,-m,\' (3) 
3 t-f \m ri - 3m n o J 

where denote the energy of the i-th pion in the rj rest frame (CM). The 
variable z lies in the interval [0 — 1], where z = corresponds to events with 
3 7T° having all the same energy while for z = 1 one 7r° is at rest and the 
remaining two are emitted back to back. 

The decay amplitude is represented at leading order in terms of a single 
quadratic slope parameter a: 

|A)oo (z) I 2 ~ 1 + 2az. (4) 

In case of pure phase space ( i.e. at leading order in ChPT) one has a = 
and the z distribution is flat from z = to z ~ 0.76 and then falls to zero at 
z = 1, see Fig. [TJ 




Figure 1: Expected z distribution for pure phase space. 



Recent measurements of a with //-mesons produced almost at rest in hadro- 
and photo-production experiments are reported in Table [TJ In the same Ta- 
ble are also shown the theoretical estimates for a previously described. The 
predicted values show differences - due to large cancellations in the ampli- 
tude computation - even quoting in some case a positive sign for a contrary 
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a 


Crystal Ball (2001) [8] 
PF.LSTTTS WAS A FQl 

WASA at COSY [10] 
Crystal Ball at MAMI-B jjjj 
Crystal Ball at MAMI-C [12] 


-0.031 ± 0.004 
n n9fi -4- n m a 

-0.027 ± 0.009 
-0.032 ± 0.003 
-0.032 ± 0.003 


ChPT / LO 
ChPT / NLO [3] 
ChPT / NLO + unit, corrections [4] 
U(3) ChPT + Bethe Salpeter [6] 
ChPT / NNLO [7] 


0.000 

0.015 
-0.014 -0.007 
-0.031 ± 0.003 
0.013 ± 0.032 



Table 1: Experimental and theoretical results for the slope parameter a. 

to the experimental evidence. A precise measurement of a therefore poses a 
significant constraint to theoretical models. We present a new measurement 
of a based on a large sample of r\ mesons produced in e + e~ collisions via the 
radiative decay (p — > 777. 



2. DA$NE and KLOE 



Data were collected with the KLOE detector at DA$NE [13J], the Frascati 
e + e~ collider, which operates at a center of mass energy W = ~ 1020 
MeV. The beams collide with a crossing angle of ir — 25 mrad, producing 
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mesons with a small transverse momentum , p^ ~ 13 MeV/c. The KLOE 
detector is inserted in a 0.52 T magnetic field and it consists of a large cylin- 
drical drift chamber (DC), surrounded by a fine sampling lead-scintillating 



fibers electromagnetic calorimeter (EMC). The DC 15| . 4 m diameter and 
3.3 m long, has full stereo geometry and operates with a gas mixture of 90% 
helium and 10% isobutane. Momentum resolution is a(p±)/p± < 0.4%. Po- 
sition resolution in r — is 150 //m and a z ~ 2 mm. Charged tracks vertices 
are reconstructed with an accuracy of ~ 3 mm. 



The EMC [16( is divided into a barrel and two endcaps, and covers 98% of 
the solid angle. It is segmented into 2440 cells of cross section ~ 4.4 x 4.4 
cm 2 in the plane perpendicular to the fibers. Each cell is read out at both 
ends by photomultiplier tubes. 
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Arrival times of particles and space positions of the energy deposits are ob- 
tained from the signals collected at the two ends; cells close in time and 
space are grouped into a calorimeter cluster. The cluster energy E is the 
sum of the cell energies, while the cluster time t and its position r are en- 
ergy weighted averages. The energy and time resolutions are respectively 



a E /E = 5.7%/ \/E (GeV) and a t = 57 ps/^/E (GeV) © 100 ps. Cluster po- 
sitions are measured with a resolution of 1.3 cm in the coordinate transverse 



to the fibers, and, by timing, of 1.2 cm/^E (GeV) in the longitudinal coor- 
dinate. 



The KLOE trigger 17] is based on the coincidence of two energy deposits 
with E > 50 MeV in the barrel and E > 150 MeV in the endcaps. Moreover, 
to reduce the trigger rate due to cosmic rays crossing the detector, events 
with a large energy release in the outermost calorimeter planes are vetoed. 

3. Event selection 

The measurement is based on an integrated luminosity of 420 pb _1 corre- 
sponding to ~ 1.4- 10 9 mesons produced. This data sample contains about 
17 millions of r\ mesons. 

The detector response to the decay of interest was studied by using the 



KLOE MonteCarlo (MC) simulation program [18[. The MC takes into ac- 
count variations in the machine operation and background conditions on a 
run-by-run basis. A MC sample for both signal and backgrounds was pro- 
duced for an integrated luminosity five times that of the collected data. In 
the MC simulation of the r\ — > 3tt° decay, the signal has been generated using 



our preliminary measurement [19[ of a = —0.027. 

We search for: — > 777 with 77 — > tt° 7r° tt° and tt° — > 77 events. To select 
the final state, we require to have seven prompt photons in the event. A 
photon is defined as an EMC cluster not associated to a DC track. We fur- 
ther require that \(t — r/c)\ < 5a t , where t is the arrival time at the EMC, 
r is the distance of the cluster from interaction point, IP, c is speed of light. 
Fig. [2] shows the photon energy spectrum. The recoil photon from the two 
body decay — ?• 777 is almost monochromatic, with E irec ~ 363 MeV and 
separated from the softer photons from tt° decay. 

All events must pass a first-level selection to filter machine background and 



an event classification procedure [18j . Events with the expected final state 
signature are selected by requiring: 
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Figure 2: Photon energy spectrum in the laboratory. (Dots: data, histogram: MC). 



• 7 and only 7 prompt photons with 21° < # 7 < 159° and > 10 MeV. 
The angle between any photon pair, 9 77 , must be > 9° to reduce split 
showers. After these selection cuts we are left with ~ 4.6 • 10 6 events. 

• A constrained kinematic fit imposing total 4— momentum conservation 
and t = r/c for each photon is performed. Input variables to the fit 
are the energies, times of flight and the coordinates of clusters in the 
EMC and the beam energies. The fit improves the photon energies 
resolution: the tt° mass resolution of ~ 15.4 MeV improves to ~ 9.6 
MeV after applying the kinematic fit. The selected events must satisfy 
the requirement P x 2 > 0.01 corresponding to x 2 < 25. After this cut we 
are left with 1.9 millions of 77 — > 37r° events corresponding to a signal 
efficiency of (40.81 ± 0.01)%. At this level, the residual background 
contamination, mainly due to K$Kl decays to neutral channels, is 
estimated by MC to be ~ 0.1%. 

• To find the best combination (among 15) of the six less energetic pho- 
tons into three 7r° a pairing procedure is applied. The procedure uses 
a pseudo-x 2 variable: 



where the invariant mass of the i photon pair, in corrispon- 

dence of the j th combination; o m is the corresponding ir° mass reso- 




j = l,2,...,15. 



(5) 
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lution parametrized, as function of the photon energy resolution: 
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(6) 



the angular resolution contribution is negligible. In Fig. [3] a data-MC 
comparison of the minimum value of the pseudo-x 2 , Xmin * s shown. 



UJ 10 ! 




Figure 3: Distribution of Xmin> use d to pair photons. (Dots: data, histogram: MC) 



The fraction of events with correctly paired photons, estimated from 
MC, is named in the following as purity, P, of the data sample. While 
we define WPf = 1 - P the wrong pairing fraction to 7r°'s. To improve 
the purity a further cut is applied: Xmin < 5- The distribution of the 
invariant mass of the two photons from 7r° decay, is shown in Fig. HI 

• After the photons pairing procedure a second kinematic fit is performed 
where the constraints on tt° and 77 mass are also imposed. For the rj 
mass we used the value 547.874±0.007 stat ±0.031 syst MeV measured 
by our experiment 20] . This fit improves the z resolution by a factor 
two. 

We define three samples with different purity applying different cuts on the 
difference of the two lowest values of x 2 , Ax 2 , as reported in Table [21 The 
resolution and efficiency as function of z are shown in Fig. [5] for the Medium 
purity sample. The reconstruction efficiency, e(z), is obtained by MC for each 
z bin, as the ratio: e(z) = N rec (z)/N gen (z) where N gen)Tec are respectively the 



generated and reconstructed events. 
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Figure 4: Invariant mass of the two photons from tt° decay after cut Xmin < 5. (Dots: 
data, histogram: MC). 
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Figure 5: Medium purity sample. Left: Resolution on the z variable. Right: Reconstruc- 
tion efficiency vs. z. 



The photon energy resolution is compared between data and MC looking 
at the distribution of AE* = E* n — E* 2 ; i.e. the difference between photons 
energy in the ir° rest frame. In Fig. |6] the distribution of AE* is plotted. 
Estimating the r.m.s. of the AE* for slices of 10 MeV in E n o a difference of ( 
1 4- 1.5)% between data and MC is observed. Consequently, the MC photon 
energies have been smeared by this amount. 

Fig.[7]shows the ratio Rae^j = (A£*)^"/(Ai?*)^. The correction improves 
the agreement between data and MC on this variable. The residual difference, 
of (0.6 ± 0.2)%, is taken into account directly in the evaluation of the result. 
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Ax 2 cut 


Samples 


Purity 


Efficiency 


N. events 


2.5 


Low 


90.4 % 


(20.07 ± 0.01)% 


948471 


5 


Medium 


95.0 % 


(12.96 ± 0.01) % 


614663 


9 


High 


97.3 % 


(7.04 ± 0.01)% 


333493 



Table 2: The three samples of different purity selected by different cuts on A\ 2 - 




A (MeV) 

Figure 6: Plot of AE* — E* — E* 2 , where E* are the 7 energies from 7r° decay in n° CM. 
(Dots: data, histogram: MC). 
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Figure 7: R A e- ratio of (AE*)f^ / {AE*)^ S vs. E n o . Dots (stars) before (after) the 
correction for the difference between Data and MC. 
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4. Measurement of the slope parameter a 

The fit to the Dalitz plot is done minimizing a log-Likelihood function 
built as follows: 

-log£(a) = - X) n i logp i (a) , (7) 

where, for each bin: are the number of reconstructed events, U{ the number 
of expected events, obtained from MC taking into account the detector reso- 
lution and WPf and weighted with l + 2az. Moreover we correct for the data- 
MC differences in the WPf. To estimate it on data we use the distribution of 
z variable reconstructed using the second best pairing combination, z x 2 . This 
distribution is fit with the superposition of the MC shapes for events with 
good and wrong pairing respectively. The uncertainty on the WPf data-MC 
difference is taken into account in evaluating the systematic error. 

The fit procedure has been tested on MC by verifying that the fit repro- 
duces in output the same input value, within the statistical error. 
To obtain the final result the fit range (0 -~- 0.7), corresponding to the region 
of the phase space in which the z distribution is flat, and the Medium purity 
sample is chosen. The fit results for the three different Purity samples are 
shown in Table O Moreover, we have applied a shift of Aa = —0.0008 on 
the slope parameter a to correct the residual data-MC discrepancy in the 
photons energy resolution, see Section [3J 





Low Purity 


Medium Purity 


High Purity 


a ■ 10 4 


-319 ±29 


-301 ±35 


-308 ± 47 




92% 


85% 


91% 



Table 3: a values from fit for different purity data samples. 

In Fig. |8] a comparison between the observed and fitted z distributions is 
shown. 

5. Systematic uncertainties 

In the following we describe the sources of systematics. For each of them, 
the fit has been repeated varying the related sources and assuming as sys- 
tematic error the difference with respect to the reference value. In Table H] 
we have summarized all the systematic errors. 



10 




Figure 8: Medium purity sample: Left: observed z distribution with the corresponding fit 
overimposed. Right: data/fit ratio as function of z. 

• Analysis cuts To control the stability of the result respect to our 
analysis cuts we have moved them independently The cut on # 77 to 
reject split showers,, was varied in the interval 6° — 18° in steps of 3°. 
The photon energy threshold was also increased from 10 MeV to 40 
MeV with a step of 5 MeV. The related systematic error is very small. 

• Energy Resolution As shown in Fig. [7J the data-MC comparison of 
(AE*) rms after correction shows a residual discrepancy of (0.6 ±0.2)%. 
While the 0.6% correction has already been applied, we estimate the 
systematics related to its uncertainty to be Aa = ±3 • 10~ 4 . 

• r) mass This systematic effect has been estimated varyin g th e rj mass 
on data by ±0.031 MeV accordingly to our measurement 20!. 

• Wrong pairing fraction For the sample used the data-MC ratio of 
WPf is 1.1 ± 0.1. As mentioned in Section H] the fit procedure takes 
into account this difference. To assign the systematic error we repeated 
the fit procedure varying the WPf within the ±10% uncertainty quoted 
above. 

• Purity As a check of the MC capability to reproduce the samples purity 
and its dependence upon z, we show in Fig. [9] the ratio between the 
number of events for the High and the Low purity sample, N High /N Low, 
as a function of z. A good agreement between data and MC throughout 
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the fit range is observed. As systematic error, we take the difference 
between the a values estimated using the Low and the High purity 
sample, see Table |3j 




Figure 9: Ratio Nhi g h/Ni ow as a function of z. (Dots: data, histogram: MC). 

• Fit range and binning The fit was repeated with different values of 
the fit range from [0-j-0.6] to [O-S-l] with a step of 0.1. This is the largest 
systematic effect. Instead we find negligible effect when changing the 
bin size by a factor 2 from 0.04 to 0.02. 



Source 


Aa 


10 4 


Analysis cuts 


-1 


+1 


Energy resolution 


-3 


+3 


r] Mass 


-2 


+6 


Wrong pairing 


-6 


+5 


Purity 


-18 


+0 


Fit range 


-29 


+20 


Total 


-35 


+22 



Table 4: Summary of the systematic errors on the slope parameter a. The total systematic 
error is the sum in quadrature of the different contributions. 
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6. Conclusions 



Using a clean sample of r\ — > 37r° decays we have measured the Dalitz Plot 
slope parameter obtaining a = —0.0301 ± 0.0035 stat toxmlsyst in agree- 
ment with other recent results of comparable precision. 
The above value is also consistent with a = —0.038 ± 0.003 stat i^oog syst 
obtained from the KLOE study of the i] — > 7r + 7r~7r° decay 21] using the the- 
oretical correlations between the two decay modes. 

Our a measurement confirms the inadequacy of simple NLO ChPT compu- 
tations and the need to take into account higher order corrections. 
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